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TIME -FIXED RENDEZVOUS BY IMPULSE FACTORING WITH 
AN INTERMEDIATE TIMING CONSTRAINT 

By Richard N. Green, Jaihes F. Kibler, and George R. Young 
Langley Research Center 

SUMMARY 

A method is presented for factoring a two -impulse orbital transfer into a three- or 
four -impulse transfer which solves the rendezvous problem and satisfies an intermediate 
timing constraint. Both the time of rendezvous and the intermediate time of' aline me nt 
are formulated as any element of a finite sequence of times. These times are integer 
multiples of a constant plus an additive constant. The rendezvous condition is an equal- 
ity constraint, whereas the intermediate alinement is an inequality constraiht. The two 
timing constraints are satisfied by factoring the impulses into collinear parts that vecto- 
rially sum to the original impulses and by varying the resultant period differences and 
the number of revolutions in each orbit. i>'Five different types of. solutions arise by con- 
sidering factoring either or both of- the two impulses into two or three parts with a limit 
of four total impulses. The impulse -factoring , technique may be applied to any two- 

impulse transfer which has distinct orbital periods. 

INTRODUCTION 

A major component of space exploration is the mission involving orbital operations. 
The maneuvers in orbit often require rendezvous. In the present study, rendezvous is 
defined as a trajectory which starts at a terminal and ends at a time-related terminal 
(ref. 1); and a terminal is a specified Cartesian position and velocity vector, or equiva- 
lently, the set of six Keplerian orbital elements. A requirement of the orbital rendezvous 
may be that the maneuver be completed at specified times (time -fixed rendezvous). The 
time -fixed rendezvous is contrasted to the trajectory which has no limitation on the ren- 
dezvous time (time-open rendezvous); Another requirement of the orbital rendezvous 
may be that the transfer trajectory satisfy an intermediate timing constraint. For exam- 
ple, to reach a particular true anomaly within a prespecified timing error is an inter- 
mediate inequality constraint on the transfer trajectory. The problem addressed in this 
analysis can thus be described as a time -fixed rendezvous with an intermediate timing 
constraint. 
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An orbital transfer between two general orbits can be achieved in two impulses. 
However, the time limitation for fixed-time rendezvous severely hinders the various ren- 
dezvous schemes which require waiting for an initial orientation before making the two- 
impulse transfer. Several authors have investigated alternative methods which require 
more than two impulses to complete the rendezvous. Van Gelder and associates (ref. 2) 
accomplished a three-impulse time -open rendezvous between coplanar orbits by splitting 
one Hohmann impulse into two equivalent impulses to obtain rendezvous with the same AV 
as the Hohmann transfer. Straly (ref. 3) investigated three -impulse time-open rendezvous 
with a circular target orbit by making a phasing impulse after arriving in the plane of the 
target orbit. Roth (ref. 4) used a bi-elliptic transfer between circular coplanar orbits to 
achieve time -open rendezvous. Bender (ref. 5) investigated three -impulse time -open ren- 
dezvous by splitting one of the impulses of the nonrendezvous optimum two-impulse trans- 
fer. Eckel (ref. 6) found a minimum-time rendezvous by a technique similar to Bender's. 
Finally, Doll (ref. 7) determined optimal multiple impulse time -fixed rendezvous trajec- 
tories, but his method requires elaborate estimates of the transfer and does not always 
yield a solution. In no instance do the rendezvous schemes consider the possibility of an 
intermediate timing constraint, and most do not consider time -fixed transfers between 
general terminals. 

The purpose of the analysis presented herein is to present a solution to the time- 
fixed rendezvous between general orbits with an intermediate timing constraint. The 
approach taken is first to determine a two -impulse orbital transfer which satisfies the 
geometry (position and velocity) of the rendezvous condition and then to adjust the timing 
in the geometry solution to satisfy the time constraints. The two -impulse geometry solu- 
tion can be any set of maneuvers which transfers from a completely specified initial orbit 
to a completely specified final orbit or rendezvous orbit. One example is the two -impulse 
transfer of McCue (ref. 8). An example of a two -impulse transfer tailored to a particular 
mission is presented in appendix A. Next, the timing in the geometry solution is adjusted 
to satisfy the intermediate timing constraint and the time of rendezvous. A coarse adjust- 
ment of the timing can be achieved by allowing the spacecraft to wait in orbit for a number 
of revolutions before proceeding to the next orbit of the two-impulse transfer. Varying 
the number of revolutions in each orbit will permit the timing to be adjusted without alter- 
ing the geometry solution. Vernier adjustments of the timing can be achieved by applying 
only part of the impulse, waiting a number of revolutions, and then applying the remainder 
of the impulse. If the two parts of the impulse are both applied at the same point and vec- 
torially summed to the original impulse, then the timing in orbit is adjusted while the 
resulting geometry is unchanged. This approach to solving the rendezvous problem is 
termed the impulse -factoring technique. It has also been called splitting, looping, phas- 
ing, and epoch changing by various authors. 
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The factoring of an impulsive velocity 
maneuver to adjust the timing is quite simple 
and is depicted in figure 1. The impulsive 
velocity maneuver is AV and is applied to 
the initial orbit to establish the transfer 
orbit. If the impulse is factored by the 
velocity factor a v and only part of the 
impulse is applied (a v AV), then an inter- 
mediate or factored orbit is established. 

This orbit has a period which is generally 
different from either the initial or the 
transfer orbit and causes a change in the time of periapsis passage. Adjustments in tim- 
ing can be made by the proper choice of the factor a v . To complete the maneuver from 
the initial orbit to the transfer orbit requires that the remaining part of the impulse, 

(l - a v )AV, be applied. Thus, an impulsive velocity maneuver can be factored to alter 
the timing while the orbital geometry is preserved. 

The full potential of impulse factoring has previously been ignored. Some authors 
have applied one factor and restricted it to lie in the range 0 to 1. This restriction 
results in no increase in the cost of the AV maneuver. However, in some cases the 
restriction does not allow a time -fixed rendezvous. In the present investigation, the 
range of rendezvous is extended by allowing factors outside the range 0 to 1. Even though 
the sum of the two resulting impulses may be larger than the original impulse, this tech- 
nique yields a time -fixed rendezvous solution that is often acceptable. In addition, this 
investigation considers using two factors to satisfy two timing constraints. One factor 
may be applied to each impulse of the two-impulse geometry solution, or both factors 
may be applied to the same impulse, whereby the impulse may be factored into three 
parts. There are three such schemes to factor a two-impulse geometry solution into a 
four-impulse rendezvous solution. In addition, there are two three -impulse rendezvous 
solutions which result from applying only one factor to either impulse. In all, five dis- 
tinct cases arise and are discussed subsequently in detail in the section "Analysis." 

An example problem using the impulse -factoring technique is presented in the sec- 
tion "Application." A further application is presented in reference 9. 

SYMBOLS 



Figure 1.- Impulse-velocity factoring. 
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e eccentricity 

f true anomaly, degrees 

H component of H 

H unit vector along angular momentum vector 

i inclination, degrees 

I,J,K,L number of revolutions relative to f° in initial, a-, transfer, and 

/3-orbit, respectively 

l',j',K',L' number of revolutions relative to f a in initial, a-, transfer, and 
j3-orbit, respectively, prior to alinement 

I a number of revolutions relative to f a prior to alinement 

I s upper bound on total number of revolutions 

m integer denoting time of alinement (see eq. (2)) 

M number of revolutions relative to f a in final orbit prior to alinement 

n integer denoting time of rendezvous (see eq. (1)) 

P orbital period, hours 

P,Q,W orbital plane coordinate system, where P is directed toward periapsis, 
Q is in orbit plane advanced from P by a right angle in direction 
of increasing true anomaly, and W completes right-handed system 

r radius from center of planet, kilometers 

? unit radius vector 

R rotation matrix from the P,Q,W to the X,Y,Z coordinate system 

t time, hours 
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time of aline ment, hours 


ta 

tr 

V 

X,Y,Z 

X',Y',Z' 

a- 

Qfy 

P 

^V 

V 
6 

6 p 

av t 

e 

c 

V 
9 
X 
A 
M 


time of rendezvous, hours 

orbital velocity, kilometers/second 

rectangular Cartesian coordinate system 

rectangular Cartesian coordinate system defined in appendix A (sketch B) 

period factor to establish a -orbit 

velocity impulse factor to establish at -orbit 

period factor to establish /3 -orbit 

velocity impulse factor to establish j3-orbit 

flight-path angle, degrees 

latitude constraint on spacecraft, degrees 

latitude of periapsis, degrees 

sum of velocity impulses, kilometers/second 

timing error at alinement, hours 

time between alinement opportunities, hours 

time between rendezvous opportunities, hours 

declination of spacecraft, degrees 

discretization variable for e 

right ascension, degrees 

gravitational constant, kilometers Vsecond^ 
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angle between two orbital planes, degrees (see appendix A, sketch C) 


Z 

r a time of first alinement opportunity, hours 

T r time of first rendezvous opportunity, hours 

( p partial sum of orbital revolutions 

\p true anomaly constraint on spacecraft, degrees 

a) argument of periapsis, degrees 

£2 longitude of ascending node, degrees 

Subscripts: 


b 

upper bound 

f 

final orbit 

min 

minimum 

0 

initial orbit 

r,h,n 

radial, heading, normal 

P,q,w 

A A A 

P,Q,W coordinate system 

t 

transfer orbit 

x,y,z 

X,Y,Z coordinate system 

x',y’,z’ 

X’,Y’,Z' coordinate system 

a 

alpha orbit 

0 

beta orbit 


i, 2, . . ., 6 first, second, . . sixth 
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Superscripts: 


a aline me nt 

i in 

o out 

Notation convention: 

closest integer to x 
increment of x 
absolute value of x 

minimum value of f(x) over the range of x 
implies 
vector x 
sign of x 


ci(x) 


Ax 



x 

sgn(x) 


ANALYSIS 

The method of impulse factoring is applied to a two -impulse orbital transfer so that 
both an intermediate timing constraint and a final rendezvous constraint are satisfied. 

The intermediate timing constraint is defined as an alinement of the spacecraft with some 
object of interest. For example, alinement could correspond to a close approach with 
another spacecraft or satellite or it could correspond to the reconnaissance of a surface 
feature. Essentially, it requires that the spacecraft reach a given true anomaly within a 
specified interval of time. In addition, the rendezvous condition requires that the space- 
craft achieve a specified position and velocity at a specified time. 

Prior to the application of the impulse -factoring method, a two-impulse orbital 
transfer which satisfies the orbital geometry of the rendezvous condition is determined. 
This "geometry" solution establishes the required position and velocity at rendezvous but 
does not establish the proper timing. The geometry solution can be determined by various 
means and is dependent upon the specific problem. Since the impulse -factoring method 
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is independent of the characteristics of the geometry solution, except that the orbital 
periods are distinct, the determination of the geometry solution is not discussed herein; 
however, an example of a typical two-impulse solution is presented in the section 
"Application.” 

The geometry solution specifies the orbital elements of the initial, transfer, and 
final orbits, which remain fixed. Quantities defined on these three orbits are denoted by 

subscripts o, t, and f, respectively. Also 
specified (fig. 2) are the true anomalies where 
the spacecraft goes "into" the orbits f* and 
where it goes "out” of the orbits f°. In addi- 
tion the true anomaly of the alinement f a 
can be determined for the specific problem 
of interest. All true anomalies represent a 
point in orbit and are defined between 0° 
and 360°. The time increment At between 
two true anomalies is denoted by its super- 
scripts. For example, AtJ° is interpreted 
as the time interval on the transfer orbit 
from fj to f° and is a positive quantity 
less than one orbital period. The actual 
Figure 2. - irue anomalies. time in an orbit is At io plus an integral 

number of orbital periods. The magnitudes of the two -impulse maneuvers are denoted 
by AV^ and AVg. These two velocity maneuvers are factored to satisfy the timing 
constraints. 

The rendezvous condition requires that the spacecraft achieve a specified position 
and velocity at a specified time. The proper position and velocity are satisfied by the 
choice of the geometry solution and the proper timing is achieved by factoring the impulse 
maneuvers. If the point of rendezvous is defined as the entry point into the final orbit t\ 
the time of rendezvous t r is the sum of the time spent in each orbit from the initial 
point f Q to the final point ff. The initial point f j, corresponds to zero time. Thus, 
the timing constraint for rendezvous is given by 

t r = r r + nrj (i) 

where n = 0, 1, 2, . . ., n b and r r and n are constants. The parameter T r is the 
time of the first rendezvous opportunity and rj is the time between rendezvous opportu- 
nities. The value of the integer variable n is bounded by n b which limits the time of 
rendezvous. The different rendezvous opportunities correspond to the times at which a 
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rendezvous is possible. For example, if the problem were to dock the spacecraft with a 
satellite for which the orbit coincided with the final orbit, then a rendezvous opportunity 
would occur each time the satellite passed through the rendezvous point ff. In this case 
r r would be the length of time for the satellite initially to achieve the true anomaly ff 
and rj would be the period of the final orbit. It might be that other mission constraints 
only permit rendezvous opportunities on alternate passages of the rendezvous point. In 
this case 77 would be twice the orbit period. 

The alinement constraint requires that the spacecraft reach a given true anomaly f a 
within a specified interval of time. If the time of alinement t a is defined as the time 
from the initial point f^ to the alinement point f a , then the inequality timing constraint 
for alinement is given by 

= ta “ ( r a + m ^) = e b (2) 

where m = 0, 1, 2, . . ., m b . The parameter T a depends upon the orbit on which aline- 
ment occurs; £ is a constant denoting the time between alinement opportunities. The 
bounds on the timing error for alinement are ± 6 ^ and the bound on the integer variable 
m is mjj. 

The orbit on which alinement occurs is not specified. However, it is required that 
alinement occur on the I a revolution. That is, f a is passed I a times before aline- 
ment occurs. It should be emphasized that I a is a specified constant and imposes a 
constraint on the solution. If I a = 0, then the alinement must occur the first time the 
spacecraft reaches f a . Although I a defines the number of revolutions prior to aline- 
ment, it does not define the alinement orbit until the number of revolutions in each orbit 
has been determined. Once this is accomplished, the alinement orbit can be determined 
by the logic of table 1 . 

The timing constraints on rendezvous and alinement are satisfied by the method of 
impulse factoring. Applying only part of the AV at f°, waiting a number of revolutions 
on the intermediate orbit, and then applying the remainder of the AV at the same place 
allows the timing in orbit to be altered while the orbital geometry is preserved. Thus, 
varying the percentage of the AV applied and the number of revolutions in the various 
orbits satisfies the timing constraints. The factored intermediate orbits are designated 
the a- and /3-orbits since they result from applying the factors a v and /3y to either or 
both AVi and AV 2 . Corresponding to each velocity factor is a period factor which 
denotes the accompanying change in period. For example, if AVj is factored by a v 
so that the a-orbit is intermediate to the initial and transfer orbits, then the resulting 
period of the a -orbit is given by 
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(3) 


B a ~ + Q! (^ > t " P 0 ) - + a AP t0 

where a is the period factor corresponding to a v . Since the constraints are expres- 
sions of time, the period factors a and /3 are more tractable than the velocity factors 
a v and 0 V . Once the period factors are determined, the corresponding velocity factors 
are easily determined (appendix B). A factor between zero and one results in two 
impulses which have the same total AV as the single impulsive maneuver. A factor 
outside this range results in two impulses which sum to a higher total AV than the 
original single impulse. If the factor is between zero and one, the transfer is said to be 
"free” (without penalty). 

The impulse -factoring method also determines the number of revolutions in each 
orbit. A revolution is defined as the number of times the spacecraft passes the maneuver 
point f° without performing the maneuver. The number of revolutions in each of the 
four orbits (initial, a, transfer, and j8) are denoted by the integer variables I, J, K, 
and L, respectively. Frequently, it is advantageous to limit the total number of revolu- 
tions prior to rendezvous. That is, 


I+J+K + Lil s 


(4) 


The two-impulse geometry solution may be factored into three or four impulses by 
considering all combinations of factoring either or both impulses into two or three parts. 
For convenience, define "bisect' 1 to denote factoring into two collinear parts that vecto- 
rially sum to the original impulse, and "trisect” to denote factoring into three collinear 
parts that vectorially sum to the original impulse. If the impulse is not factored, it is 
said to be applied in "full." Limiting the total number of impulses to four or less allows 
five different types of solution: (1) bisect-full (BF), (2) full-bisect (FB), (3) bisect -bisect 

(BB), (4) trisect-full (TF), and (5) full-trisect (FT). The first two types are three-impulse 
solutions and the remainder are four -impulse solutions. It has been assumed that the 
full -full or geometry solution does not satisfy the time constraints. 

The rendezvous problem with an intermediate time constraint as proposed herein 
does not have a unique solution. In fact, it is possible for all five types of solutions to 
satisfy the time constraints. The desired solution is the one which minimizes the sum 
of the impulses. For example, in a solution of the bisect -bisect type, the sum of the 
impulses is given by 


AV-ji — 


a. 


AVi + 


T - a, 


AVi 


K AVo + 


1 -Pv 


AVr 


(5) 


Thus, the problem consists of determi ni no* tVio cof r%f infafTAr Xro vioVdoo T T IT T 

' - " ~ — **'**C> “‘ 1V/ WV V* iUMll/U Xj O y i-l y 1X1 

n and the factors ot and which satisfy the rendezvous and alinement constraints at 


10 



the least cost for the given geometry solution. To determine this solution, the variables 
I, J, K, and L are systematically varied within the following limits: 

I ^ Wn = 0 

J = Jmin > 0 

K ? K min a 0 

L = L min > 0 
I+J + K + L = Ig 

where the minimum values of I, J, K, and L are imposed on the solution. The min- 
imum values of J and L must be greater than zero to avoid eliminating their respec- 
tive orbits. Next, for each set of I, J, K, and L, the best values of a, fi, m, and n 
are determined for each of the five different types of solutions. In this manner, the type 
of solution and the values of I, J, K, L, m, n, a, and /3 which minimize the sum 
of the impulses are determined. 

In the remainder of this section the necessary equations for each of the five types 
of solution are developed. 


Bisect -Full Solution 

The bisect-full type of solution is a three-impulse solution where AV^ is factored 
by ct v to establish the a -orbit, which follows the initial orbit and precedes the transfer 
orbit. The AV 2 impulse is unaltered and applied in ''full”; that is, the spacecraft tra- 
verses first the initial orbit, then the a -orbit, the transfer orbit, and the final orbit. The 
rendezvous time is satisfied by the a -factor which also dictates the alinement time. If 
the alinement time is within bounds, then the solution is acceptable, which means that the 
solution satisfies the two time constraints. This result is not intended to infer, however, 
that the solution is complete since other acceptable solutions may yield a lower total AV. 
If the alinement time is not within bounds, then the bisect-full solution is unacceptable. 

For the bisect -full type of solution, the spacecraft proceeds from f^ on the initial 
orbit to the true anomaly of exit f°. It then makes I complete revolutions in the initial 
orbit relative to f° before transferring to the a -orbit. Since the transfer in and out of 
the a -orbit occurs at the same point ^ = fg), the spacecraft makes J complete revolu- 
tions in the a -orbit before transferring into the transfer orbit. It then traverses to f°, 
makes K complete revolutions relative to f°, and transfers to the final orbit at f| 
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which is the rendezvous point. Thus, equating the rendezvous time (eq. (1)) and the time 
in orbit until the spacecraft arrives at f^ yields 

Atj, 0 + IP 0 + JP a + Atj° + KP t = T r + nri (6) 

where n = 0, 1, 2, . . n b . Substituting equation (3) for Pq, into equation (6) yields 

Atj, 0 + IP 0 + J(P 0 + a AP to ) + Atj° + KP t = r r + n V 

For a given value of I, J, and K, 

Cja + C 2 = n7? (7) 

where 


Cl = J AP t0 


C 2 = (I + J)P G + KP t + At™ + At{° - r r 


and the C's are constants. Thus, the value of a is given by 


n7j - Co 
( X = -= — — 

Cl 




and the total cost is 


AVrp = 


a. 


AV 1 


1 - a. 


AVl + AV r 


(9) 


where a? v is determined from a (appendix B). 

Next, the defining equation for the alinement is written. Since this equation takes 
a different form for alinement on different orbits, each possibility is considered in turn. 
For a given value of I, J, and K, the alinement orbit is determined by table 1. 

(a) Initial orbit . - Assume alinement on the initial orbit. Equating the time of aline- 
ment to the time of arrival of the spacecraft at the alinement point yields 


Atjf- + i r p Q + e = r ao + m£ (10) 

where m = 0, 1, 2, . . ., m b and Atjf is the time increment on the initial orbit from 
fj, to fg. The time error in the spacecraft arrival at the alinement point is denoted 
by e. Solving equation (10) for e and minimizing over m gives 
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e - min 
m 


mC + C 3 


(m = 0, 1, 2, . . ., m b ) (11) 


where 


c 3 - T a,o - < a - I,p 0 

The parameter I' is the number of revolutions referenced to f a before the spacecraft 
reaches the alinement point. Summing the angular traverse to the alinement point and 
equating it to the required traverse yields 

Af* a + 360° I a = Afj> a + 360° I' 


or 


I’ -k 

Now, if e = Sjj, the solution is acceptable and AV T can be minimized over n, where 
n = 0, 1, 2, . . n b . 

(b) a -orbit. - Assume alinement on the a -orbit. The defining equation for alinement 
is 

Atj>° + IP D + Atjf + J'Pq, + e = T a>a + (12) 

The values of Atjf and r a a require a knowledge of the a -orbit which is yet to be 
determined. To overcome this difficulty, assume a linear transition of these quantities 
with respect to a between the initial and transfer orbits. Thus, the following approxi- 
mations are made: 


* At° a + aLtl* - At° a ) 

(13) 

* T a,o + “( T a,t " T a,o) 

(14) 


Substituting equations (3), (13), and (14) into equation (12) yields 

At* 0 + IP 0 + At° a + o^Atf - At° a ) + J'(P 0 + a AP t0 ) + e = r ao + o(r a>t - + m? 

or 


Cga; + C4 + e = m£ 


(15) 
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where 


c 3 = J ' AP to * ‘“I* ' A C - ’■a.t + T a,o 

C4 = (I + J )P 0 + At Q + At Q - r a 0 

If the position of alinement is assumed to be relatively insensitive to changes in orbit so 

that Af^ a + 360° I a is a good approximation to the angular traverse prior to alinement, 

then 

Afj, a + 360° I a ~ Af a ° + 360° I + Af jf + 360° J' 

In 03 r 

Assuming that ~ Af 0 , which is consistent with the previous assumption, and solving 

for J' yields 


J’ 




A 

+ ^o 
360° 



where the operator "ci" denotes closest integer. Solving equation (15) for e and elimi- 
nating a with equation (8) gives 


e = min 
m,n 




(16) 


where n = 0, 1, 2, . . ., n^ and m = 0, 1, 2, . . ., If e § e^, then the solution is 

acceptable. Usually only one set of m and n will satisfy e = e^. However, if more 
than one acceptable solution exists, one can determine the acceptable values of m and n 
which yield the least cost according to equation (9). 

(c) Transfer orbit. - Assume alinement on the transfer orbit. The defining equation 
for alinement is 

At Q ° + IP Q + JPqt + At| a + K P^ + e = T a ^ 

Substituting for P a and collecting terms yields an equation of the form 

Cgtj + C4 + 6 = m£ (17) 
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where 


C 3 = J *P to 

C 4 = (I + J)P 0 + K’p t + At^° + At{ a - T a> t 


K' = ci 


I a 


I - J - 


Af n° + Af t ia - Af^ 
V 360° / 


The simultaneous solution of equations (8) and (17) is given by equation (16). 

(d) Final orbit. - Assume alinement on the final orbit. To this point it has been 
accepted that the rendezvous condition is satisfied at the entry point into the final orbit. 
However, once the rendezvous condition is satisfied, it is continuously satisfied unless 
the orbit is altered. Since the final orbit remains unaltered, alinement on the final orbit 
is allowed to occur on any successive passage of f a provided it is on the I a th revolu- 
tion and the time of alinement does not exceed t a = r a f + m b £ + e b . If the number of 
revolutions relative to f a in the final orbit prior to alinement is denoted by M, the 
defining equation for alinement is then 

Atjp + IP Q + JP a + AtJ° + KP t + At| a + MP f + e = + m? 

Substituting for P a and collecting terms yields an equation of the form 

Cgor + C 4 + e = m£ (18) 


where 


C 3 = J AP t0 

C 4 = (I + J)P Q + KP t + MP f + Atjp + At\° + At| a - r a<f 


M = ci 


I - J - K 


Af, 


lO 


Af t 


+ Af| a - 


Af 


ia^ 


360 c 


The simultaneous solution of equations (8) and (18) is given by equation (16). 

For the bisect -full type of solution the alinement can occur on four different orbits. 
Each of these possibilities has been fully investigated and results are summarized in 
table 2. 
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Full-Bisect Solution 


The full -bisect type of solution is a three -impulse solution where AV 2 is factored 
by /Sy so that the /3-orbit follows the transfer orbit and precedes the final orbit. The 
AVj impulse is unaltered. This type of solution is similar to the bisect -full type with 
the difference being the factoring of AV 2 as opposed to the factoring of AVj. Since 
the derivations of the governing equations are so similar, they are not developed here but 
are summarized in table 3. 


Bisect-Bisect Solution 

The bisect -bisect type of solution is a four -impulse solution where AVj is factored 
by a v and AV 2 is factored by ^ in such manner that the order of the orbits is as 
follows: initial orbit, a -orbit, transfer orbit, /3-orbit, final orbit. In this type of solution 
the two factors (o' and /3) are so determined that the two time constraints are satisfied. 

The solution is determined by the simultaneous solution of the equations governing 
the alinement and rendezvous conditions. The rendezvous condition is 

Atjp + IP 0 + JPq, + At{° + KP t + LiP^g = r r + n?/ 

An expression for the period of the /3-orbit is 

P/3 = p t + 0 AP ft ( 1Q ) 

which is similar to the expression for P a (eq. (3)). Substituting these expressions for 
the periods in the rendezvous equation yields 

CjO! + C 2 /3 + C 3 = nr? (20) 


where 


Cl - J AP t0 


C 2 - L APjj. 

C 3 = (I + J)P G + (K + L)P t + Atj? + Atj° - T r 
The cost of the solution is given by 


AVt — 


a. 


AVi 


1 - a. 


AVi 




A Vo + 


1 -/ 3 , 


AVfl 


Next, the governing equations for alinement are derived. 


( 21 ) 
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(a) Initial orbit .- Assume alinement on the initial orbit, The alinement condition is 
given by equation (11) as 


e = min 
m 


m£ + C4 



where 


c 4 - T a,o " " W> 

If e = e^, then the solution is acceptable and AV-p is minimized subject to equation (20). 
From equation (20), 


Q! 


mj - C 2 /3 - C 3 


Cl 


( 


and substituting into equation (21) with the assumption that a v ~ a and j3y ~ 0 gives 


AV-p ~ 


n V - C 2 0 - C 3 


AVi + 


Cj - nrj + C 2 $ + C 3 


Cl 


AVi + 


av 2 + 


1 - ^ av 2 


It is desired to minimize AV-p over n and 0. For a given value of n, AV-p 
is a piecewise linear function of 0 with corners (that is, discontinuous first derivatives) 
at the points where each of the four terms equal zero. Since AV-p is a piecewise linear 
function, its minimum must occur at one of the four corners created by the absolute- 
value operator. Another possibility exists if there is a range of 0 in the closed inter- 
val 0 to 1 so that a is also in the same range. In this case, the factoring is "free” and 
AV-p = AV^ + AV 2 over the closed interval 0 to 1. Thus, there exists a family of solu- 
tions which minimize AV T . It can be shown that one of these solutions corresponds to 
one of the solutions obtained by setting each term of AV-p equal to zero. It follows that 
for either case AV-p is minimized by either a = 0, a = 1, /3 = 0, or 0=1 and the 

four -impulse bisect -bisect solution degenerates to one of the types of three-impulse solu- 
tions. All four cases must be investigated by using equation (20) to determine the other 
factor; and AV-p is minimized over n = 0, 1, 2, . . ., n^. Since a v = a and 0y = 0 
at the two points zero and one, the approximation is exact for the four degenerate cases. 

(b) Q! -orbit . - Assume alinement on the a -orbit. The alinement condition is given 
by equation (15) as 


C40! + Cg + e = m£ 


( 22 ) 
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where 


C4 - J AP^ 0 + Atj - At Q r a ^ + r ao 
C 6 = (I+J’)Po +At o 0 + A C- T a,o 



The solution is obtained by minimizing AVVp subject to equations (20) and (22) where 
-e^ = e = e b , n = 0, 1, 2, . . n b , and m = 0, 1, 2, . . m^. The solution is simplified 

by discretizing e as e = where X = -1, 0, 1. Thus, solving equation (22) for a 
and equation (20) for /3 yields 

AV t = min f lad AVj + |l - Qf v | AVj + 1/sJ AV 2 + 1 - ^ AVjjl (23) 

n,m,xL 1 1 1 1 J 


where 

(“ - 0 ! V ) 

“ ^v) 

n - 0, 1, 2, . . n b 
m = 0, 1, 2, . . m b 
X = -1, 0, 1 

and a v and /3y are derived in appendix B. 

(c) Transfer orbit .- Assume alinement on the transfer oi'bit. The alinement is 
given by equation (17) as 

C^a + Cq + e = m£ (24) 
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where 


C 4 - J AP to 

C 6 = (I + J)P 0 + K ’ P t + K + At t a - T a,t 
K' = ci I a - I - J - 

The solution to equations (20) and (24) is given by equation (23). 

(d) /3-orbit. - Assume alinement on the /3-orbit. The defining equation for alinement 
is 

At^ 0 + IPq + JPq* At|^ + KPj- + At^ a + L + 6 = ^a,/3 

where 

At]f ~ At£ a + /3 (At - At£ a ) 

T a,/3 ~ T a,t + ^( T a,f ~ T a,t) 

Making these substitutions along with equations (3) and (19) yields 

C A a + C 5 /3 + C 6 + e = (25) 

where 

C 4 = J 4P to 

C 5 = Atf - it? 1 + L’ AP ft - T a>f + r a t 

C 6 3 (I + J)P 0 + (K + L')P t + 4t‘° + 4tj° + 4tt“ - T a _, 

Afp 0 + Af| 
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In general, equations (20) and (25) can be solved simultaneously for a and 0. Thus, 
the solution is 


AVrp - min J 
n,m,X I 


a. 


AVi 


1 - ajAVi + 


AVo + 1-0, 


AV, 


where 


‘1 

- C 3J 

l C 5 

+ 1 

(mC - < 

“6 ” Xe b) 

lC 2 

c l< 

fC. 

\ 

2 ’ C 5) 

l 


0 = 




(a - <x v ) 


( 26 ) 


(27) 


(28) 


n — 0 , 1 , 2 , . . n b 
m = 0 , 1 , 2 , . . m b 
X = -1, 0, 1 

and the relation C 4 = Cj has been incorporated. 

This solution is well behaved unless = 0 or C 2 = C 5 . The constant 
Cj = J AP |- 0 must be nonzero since J > 0 and AP to # 0 because the geometry solu- 
tion is restricted to orbits with distinct periods. If C 2 = C 5 , then from equations (20) 
and (25) the minimum value of e is 

e = minjj^m? - C 3 ) - (m£ - C 6 )|j (29) 

where n = 0 , 1 , 2 , . . ., n b and m = 0 , 1 , 2 , . . ., m b , and 


a _ ( n? 7 - c 3) - C 2 ^ 

C 1 


( a - a v) 


(30) 


If € § 6 j y then the solution is acceptable and AV-p is minimized. From equation (30) 
AV«p can be expressed in terms of four continuous linear functions of either a or 0. 
Since the minimum of AV-p must occur at one of the corners created by the absolute 
operator, either a = 0, a = 1, /3 = 0, or 0=1. The minimum of AVp> is chosen from 

these four possibilities. 
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(e) Final orbit. - Assume alinement on the final orbit. The defining condition for 
alinement is given by 

Atj, 0 + IP Q + JPq, + AtJ 0 + KP t + LP^ + At} a + MP f + e = r a?f + m £ 

Substituting for Pq, and P^ and collecting terms yields 

CjO! + C 2 0 + C4 + € = m ^ (31) 


where 


C4 = (I + J)P G + (K + L)P t + MP f + Atj 5 + Atf° + Atf - T^ f 


M = ci 


Is -I-J-K-L- 


' Af^° + AfJ° + Af| a - Af Jf 
V 360® j 


From equations (20) and (31) the minimum value of £ is 

e = mini (nv7 - C 3 ) - (m£ - C4) 1 
m,nl_ J 


where n = 0, 1, 2, . . ., n b and m = 0, 1, 2, . . ., m b . If the solution is accept- 

able and AV t is minimized subject to equation (20); that is, 


nr} - C 2 /3 - C 3 
Cl 


(« ~ «v) 


(32) 


From equations (32) and (26), AV T can be expressed as the sum of four continuous 
linear functions of either a or 0. Since the minimum of AV T must occur at one of 
the corners created by the absolute operator, either a = 0, a ~ 1, 0 = 0, or 0=1. The 
minimum of AV-p is chosen from these four possibilities. 

For the bisect-bisect type solution the alinement can occur on any of the five orbits. 
Each case has been investigated and is summarized in table 4. 

Trisect -Full Solution 

The trisect-full type of solution is a four -impulse solution in which AV^ is factored 
into three parts by a and 0 and AV 2 is applied in full. The order of the orbits are 
the initial orbit, a -orbit, 0-orbit, transfer orbit, and final orbit. 
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The period of the a -orbit is the same as previously defined; that is, 


P(* - Po + a (Pf Po) - Po + a AP to 
and the period of the /3-orbit is 

P/3 “ P« + ^(Pf " Pa) 
or 

P^= P Q + (« + /3 - a/3)AP to (34) 

where P a was eliminated by use of equation (33). Thus, the a -factor bisects the period 
difference between the initial and transfer orbit, and the /3-factor bisects the period dif- 
ference between the a- and transfer orbits. The result is to trisect the AVj impulse. 
The defining equation for rendezvous is 

At™ + IP 0 + JP a + LPp + At™ + KP t = r r + nr] 

Substituting for P a and P^ (eqs. (33) and (34)) and collecting terms yields 

Cj^af + C 2 /3(l - a) + C 3 = nrj (35) 


where 


C 1 = (J + L)AP t0 
C 2 = L AP to 

C 3 = (I + J + L)P 0 +■ KP t + Ato° + At™ - r r 


Equation (35) can be solved for either a or /3 in terms of the other factor; that 
is, 


a = 


n - Co/3 - Co 

“Cf - 


(“ - a v) 


(36) 


or 


n?j - Cja - C 3 

C 2 (l - a) 


(e-fiv) 


(37) 
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Q 

From equation (36) note that a value of 0 = = J + L cannot be tolerated. The reason 

for this restriction becomes clear by examination of the time spent in the a - and 0-orbits; 
that is, 


JP a + ^P^ - JPq> + L P<2 + /^(Pt “ P<x) 


= (J + L)p t 


Thus, the time spent in the a - and 0-orbits is independent of a if 0 = - ~ — - and the 

-Li 

rendezvous condition cannot be satisfied. Essentially, the problem is that the a -orbit 
causes a phase shift and the 0-orbit causes an equal and opposite shift relative to the 
period of the transfer orbit. From equation (37) note that a = 1 cannot be tolerated, 
Physically, a value of a = 1 results in the a -orbit coinciding with the transfer orbit 
and there is no period difference for 0 to factor. For the acceptable solutions the cost 
is given by 


AV-p = cu v AVj + 0 V 


O', 


AVi + 1 - 0 v j 1 - a v AV X +AV 2 


(38) 


Next, the governing equations for alinement are derived. 

(a) Initial orbit . - Assume alinement on the initial orbit. The alinement condition 
is given by equation (11) as 


€ = min 
m 


m£ + C4 


(m = 0, 1, 2, . . ., m b ) 


where 


c 4 ~ T a,o At o " I a p o 

If «?e b , the solution is acceptable and AV T is minimized subject to equation (35). As 
before AVp can be expressed as the sum of three continuous functions of either a 
or 0. Since the minimum of AVp must occur at one of the corners created by the 
absolute operator, either a = 0, 0 = 0, or 0=1. The minimum of AV-j is chosen 

from these three possibilities. 

(b) 01 -orbit . - Assume alinement on the a -orbit. The alinement condition is given 
by equation (22) as 

C4Q! + C5 + e = m£ (39) 


where 


C4 - J AP bo + At^ - Atp 3, - r a t + T a Q 
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C 5 = (I+J')P 0 +Ati° + < a -r aj0 


J' = ci 


Ia-I- 


Af? + Af° a - AfJ? 
V 360° / 


Solving equation (39) for a gives 


‘ Xe b " C *> 

a = 7 ^ * 

c 4 


(a - a v ) 


(40) 


Thus, the solution is obtained by minimizing AV T where a is given by equation (40); 

0 is S iven b y equation (37); and n = 0, 1, 2, . . n b ; m = 0, 1, 2, . . m b ; \ = - 1 , 0, 1 . 

(e) ft-orbit .- Assume alinement on the j3-orbit. The alinement equation is given by 

At? + IP 0 + JP a + At? + L’P^g + e = T a ^ + m£ (41) 

The values of At^ and r a ^ require a knowledge of the a - and /3-orbits which are yet 
to be determined. This difficulty can be overcome by making the following approximations: 

At? * At£ a + (a + 0 - <*/3)(At? - At° a ) 


T a,/3 * T a,o + (« + 0 - «0 )(t m - r aj0 ) 

Substituting these approximations, along with equations (33) and ( 34 ), into equation (41) 
yields 


where 


C 4 a + Cg|S(l - a) + C 0 + e = m£ 


C 4 - (J + L’) AP t0 + At? - At£ a - V + r a o 
C 5 = L* AP t0 + At? - At° a - T a t + r a>0 
C 6 = (I + J+L')P 0 + At? + At° a -T a>0 


L’ = ci 


la - I - J - 


Af to + M 

360° 




(42) 
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(« ~ a v) 


(43) 


Solving equations (35) and (42) for a yields 


a = 


_ (n V - C 3 )C 5 - (mg - C 6 )C 2 + C 2 Xe b 
c l c 5 " C 2 C 4 


Thus, the solution is obtained by minimizing AV T where or is given by equation (43); 

$ is given by equation (37); and n = 0, 1, 2, . . n b ; m = 0, 1, 2, . . m b ; X = -1, 0, 1. 

(d) Transfer orbit. - Assume alinement on the transfer orbit. The alinement equation 
is given by 

Atjp + IP 0 + JP a + LP^ + Atj a + K'P t + e = T a t + mg 


Substituting for P a and P^ and collecting terms yields 
Cjdi + C 2 / 3(1 - oi) + C 4 + c = mg 


(44) 


where 


C 4 = (I + J + L)P 0 + K’P t + AtJ, 0 + At{ a - r a t 


K' = ci 


I - J - L - 


Afp° + Af{ a - Afjf 
v 360° ) 


and Cj and C2 are given by equation (35). It is readily seen that equations (35) 
and (44) cannot be solved simultaneously for a and / 3 . Thus, solving these equations 
for e and minimizing over m and n gives 

e = min/ -(n^ - C3) + (mg - C 4 
m,nlj v x 

where n = 0, 1, 2, . . ., n^ and m = 0, 1, 2, . . ., m^. If e £ the solution is accept- 
able and AV-y is minimized over the factors subject to equation (35). Here, again, the 
minimization of AVVp can be reduced to three possibilities, a = 0, /3 = 0, and 0 = 1 . 

(e) Final orbit .- Assume alinement on the final orbit. This solution is very similar 
to the solution for alinement on the transfer orbit. Therefore, it is not repeated but is 
summarized in table 5 along with the five cases of the trisect -full type of solution. 

Full -Trisect Solution 

The full-trisect type of solution is a four-impulse solution where AV^ is applied 
in full and AV2 is factored into three parts by a and 0. The order of the orbits is 
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the initial, transfer, a-, fi-, and final orbit. This type of solution is similar to the trisect - 
full solution with the difference being the trisecting of AV2 as compared with trisecting 
AVj. Since the derivation of the governing equations are so similar, they are not devel- 
oped here but are summarized in table 6. 

APPLICATION 

As an example of the impulse -factoring technique, consider a satellite mission about 
the planet Mars. Assume the satellite is initially in a known orbit. A sequence of impul- 
sive maneuvers is desired that will cause the spacecraft to fly over a selected surface 
feature for reconnaissance and then establish a synchronous orbit in such a way that the 
surface feature is directly beneath the spacecraft at each periapsis passage. This type 
of mission might involve surface mapping or it might involve establishing a proper orbit 
from which to launch a surface probe. 

The geometry of the mission is presented in figure 3. The Keplerian orbital ele- 
ments of the initial orbit and the right ascension of the surface feature are known at some 
epoch. The first task is to define a two-impulse geometry solution that establishes the 
desired final orbit. The conditions imposed on the final orbit are that it have a synchro- 
nous period and that periapsis be directly over the latitude of the surface feature. In 



Figure 3>- Mission geometry. 
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addition, the altitude of periapsis is fixed* Thus, the synchronous period and the periapsis 
altitude define the shape of the orbit, that is, the semimajor axis a and the eccentric- 
ity e. The condition of periapsis over the correct latitude leads to a relationship between 
the orientation angles i, £2. Since the final orbit is not fully defined, the geometry 
solution is not unique and various two-impulse solutions are possible. A specific geom- 
etry solution can be defined, however, by choosing the magnitude and direction of the first 
impulse, the position in the initial orbit at which it is applied, and the true anomaly of 
entry into the final orbit. These five independent quantities along with the quantities 
already specified uniquely define the geometry solution. (See appendix A.) Since the 
choice of the five independent quantities is arbitrary, these quantities can be varied to 
minimize the sum of the impulses for the geometry solution. (See ref. 10.) This mini- 
mization gives the following solution: 



Initial orbit 

Transfer orbit 

Final orbit 

a 

*20762 

21242 

*20428 

e 

*0.77524 

0.76285 

*0.76045 

i . . . . ; 

*33.20 

33.86 

33.24 

CO 

*34.38 

30.32 

28.17 

a 

*104.85 

107.39 

107.93 

p 

*25.23 

26.11 

*24.62 

fi 

*0.00 

213.51 

125.87 

f° 

211.58 

124.18 


f a 

*-6.17 

-2.63 

*0.00 

Af io 

211.58 

270.67 


Af ia 

*353.83 

143.86 

234.13 

Af oa 

142.25 

233.19 


At i0 

21.04 

6.30 


At ia 

*25.20 

4.22 

22.55 

At oa 

4.16 

24.03 


T a 

*7.86 

7.99 

8.06 

6 P * * * 

*18.01 

16.33 

*15.00 


*Fixed quantity. 


The true anomaly of aline ment f a depends on the orbit. For this example it is defined 
as the true anomaly directly above the latitude of the surface feature which is 15°. Like- 
wise, the time constant for alinement r a depends on the orbit and is defined as the time 
required for the surface feature to rotate beneath the orbit, that is, to rotate to the point 
of alinement. Since the object of alinement is a surface feature, the time between aline- 
ment opportunities 5 is the rotational period of Mars (24.62 hours) which is also the 
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time between rendezvous opportunities 77 since the final orbit is synchronous. The time 
constant for rendezvous r r causes the surface feature to be directly beneath the space- 
craft at periapsis passage on the final orbit. Numerically, it is equal to the time required 
for the surface point to rotate beneath the final orbit minus the time required for the 
spacecraft to pass from ff to periapsis (recall that rendezvous occurs at f|), that is 

T r = T a,f " At f a 
r r = 8.06 - 22.55 

r r = -14.49 or 10.13 

Since rendezvous cannot occur at a negative time, and since the time of rendezvous is 
t r = r r + n?7, where n = 0, 1, 2, . . n b , the time constant is modulo 77 or r r = 10.13. 
As a matter of interest, the latitude of periapsis 6p on the three orbits is shown. The 
only other quantities needed for the impulse-factoring technique are: 


AVi = 0.03467 

Lmin = 

AV 2 = 0.02688 

I a = 4 

% = 0,5 

I s = 8 

Imin ~~ 1 

mjj = 12 

J min = 2 

n b =ll 

Kmin = 1 



The sum of velocity impulses is minimized for each type of solution and the results are 
presented below: 


Type of 
solution 

I 

J 

K 

L 

m 

n 

a 

0 

av t , 

km/sec 

Bisect -full 

1 

2 

1 

— 

5 

5 

2.336 


0.14369 

Full -bisect 

1 

— 

1 

2 

5 

5 

— 

-0.790 

.09995 

Bisect -bisect 

1 

3 

1 

3 

5 

9 

1.676 

1.206 

.11635 

Trisect -full 

1 

3 

1 

3 

5 

9 

1.623 

4.188 

.24715 

Full -trisect 

1 

2 

1 

4 

5 

9 

-.586 

.936 

.09036 
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None of the five types of solutions are free since all have factors outside the range 0 to 1. 
The least expensive solution, however, is the "full -trisect” type of solution where aline - 
ment occurs on the a -orbit (note f a is passed twice on the transfer orbit) after approxi- 
mately 5 Martian days and rendezvous on the ninth day. In order to satisfy the alinement 
constraint, a is required to be negative and results in additional cost relative to the 
cost of the geometry solution. However, this is the least expensive way to satisfy the 
rendezvous condition with an intermediate timing constraint with the given geometry solu- 
tion. The flow diagram for this example is given in figure 4(a). 

In the preceding numerical example the cost of the geometry solution (aVi + AV 2 ) 
was minimized by varying the five independent quantities associated with the solution and 
then the impulse -factoring technique was applied. If the resulting factors are outside the 
range 0 to 1, then the total cost AVt will exceed the cost of the geometry solution. For 
this situation the total cost may be reduced at the expense of additional calculations by 
minimizing AV-p over the five independent quantities. This approach involves applying 
the impulse -factoring technique at each step during the minimization process instead of 
only once as shown in figure 4(b). The result is to change the geometry solution in such 
a way that the factoring technique is less expensive. The geometry solution resulting 
from minimizing AV-p is given in the following tabulation: 



Initial orbit 

Transfer orbit 

Final orbit 

a 

*20762 

21563 

*20428 

e 

*0.77524 

0.76868 

*0.76045 

i 

*33.20 

33.34 

33.04 

U) 

*34.38 

30.33 

28.34 

a 

*104.85 

106.54 

106.98 

p 

*25.23 

26.70 

*24.62 

f l 

*0.00 

230.14 

113.58 

f° 

227.51 

111.95 


f a 

*-6.17 

-2.24 

*0.00 

Af io . . 

227.51 

241.81 


Af ia 

*353.83 

127.62 

246.42 

Af oa 

126.32 

245.81 


At i0 

22.81 

3.89 


At ia 

*25.20 

2.42 

23.15 

At oa 

2.38 

25.23 


T a 

*7.86 

7.96 

8.01 

5p 

*18.01 

16.12 

*15.00 


"“Fixed quantity. 
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(a) Minimum geometry cost solution. 


(b) Minimum total cost solution. 


Figure 4.- Flow diagrams of impulse- factoring solution. 












The cost of the geometry solution is AV^ = 0.03906 and AV 2 = 0.02719 for a total cost 
of 0.06625 km/sec. This result is slightly greater than that of the first example which 
cost 0.06155 km/sec. However, in the second example the geometry solution has been 
changed to reflect the factoring process. The results of the five types of solutions are 
presented below: 


Type of 
solution 

I 

J 

K 

L 

m 

n 

a 

0 

av t , 

km/sec 

Bisect- full 

1 

2 

1 

— 

5 

5 

1.186 


0.07969 

Full- bisect 

1 

— 

3 

2 

5 

7 

— 

0.868 

.06625 

Bisect-bisect 

1 

2 

1 

2 

5 

7 

1.000 

.868 

.06625 

Trisect-full 

1 

3 

1 

3 

5 

9 

.973 

-63.013 

.20688 

Full-trisect 

1 

2 

1 

2 

5 

7 

.000 

.868 

.06625 


It is seen that the "full-bisect," "bisect-bisect, " and the "full-trisect" solutions are iden- 
tical and that they are free. Thus, selecting a geometry solution which minimizes AV^ 
instead of (AVj + AV 2 ) reduces the total cost from 0.09036 to 0.06625 km/sec. The 
reduction in cost is not always this appreciable; however, when the cost of the factored 
solution greatly exceeds the cost of the geometry solution, the potential for a cost reduc- 
tion exists by minimizing AV-j'. 

CONCLUDING REMARKS 

A method has been presented for factoring a two-impulse orbital transfer into a 
three- or four-impulse transfer which solves the rendezvous problem and satisfies an 
intermediate alinement constraint. Five types of solutions exist and depend upon the fac- 
toring scheme employed. The equations governing each solution have been derived. 

The impulse -factoring technique has many advantages which make it desirable. The 
rendezvous problem is simplified by first solving the orbital geometry transfer and then 
satisfying the time constraints by factoring the velocity-impulse maneuvers. The method 
generates a number of different solutions which satisfy the rendezvous condition and the 
alinement constraint. This set of solutions is finite for a given geometry transfer, and 
the minimum velocity solution can be chosen from the set. The method is easily pro- 
gramed and circumvents many of the mathematical and computational problems associated 
with more classical approaches. 

An application of the impulse- factoring technique has been presented. An example 
problem was solved by minimizing the sum of the two impulsive velocity maneuvers of 
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the geometry transfer and then factoring the impulses to satisfy the time constraints. 
This rendezvous solution required more velocity increment AV than the geometry 
transfer. Therefore, the example problem was reworked to minimize the AV of the 
rendezvous solutions by varying the geometry transfer, resulting in a reduction of AV. 
It is concluded that if the rendezvous solution costs more than the geometry transfer, it 
is better to determine the geometry transfer which reflects the factoring process than 
merely to factor the geometry transfer of least AV. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., April 4, 1974. 
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APPENDIX A 


TWO-IMPULSE ORBITAL TRANSFER TO ESTABLISH a, e, 6 

The necessary equations for defining a two-impulse transfer from a known Keplerian 
initial orbit to a partially specified final orbit are derived in this appendix. The initial 
orbit is defined by a 0 , e 0 , i 0 , cu 0 , fl 0 ; and the final orbit is partially specified by a^, e^, 
and 6 where 6 is the latitude of the spacecraft when the true anomaly has the value \f/. 
The two velocity maneuvers are defined as (see sketch A): 

AV7= (AVl ?r , AVi )h , AVi >n ) 

(AV 2>r , AV 2> AV 2?n ) 


w 



The problem to be solved can be stated as follows: 

Given: a 0 , e Q> i Q , o> 0 , 0 0 , aj, e f , 5, fi, fg, AV l r , AV l h , AV 1>n , f f 
Find: AV^ and the complete description of the transfer and final orbits. 

The parameter fg is the true anomaly in the initial orbit where AV j is applied 
and f| is the true anomaly in the final orbit corresponding to the AVlj maneuver. 
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APPENDIX A - Continued 

Although the parameters fg, AV^ r , AV^ h, AV^ n , fj could be considered independent 
variables and AVj + AV 2 could be minimized, this problem is not addressed here. 

A ^ A 

The first transfer is straightforward. Consider the initial orbit in the P,Q,W 
coordinate system (sketch A). The magnitude of the first maneuver is simply 

iV l - ( aV l,r + aV l,h + * V 1,J ^ 

and the elements of the resulting transfer orbit can be established by computing its com- 
ponents of position (r) and velocity (V) in the P,Q,W system, rotating these components 
to the original X,Y, Z system, and finally converting them to orbital elements. These 
calculations are as follows: 


*o( x - e o) 

1 + e 0 cos fg 



sin y Q = 


e o sin fg 


(l + 2e 0 cos fg + e^) 1 


cos y Q = [l - sin 2 

r p = r 0 cos f 0 

r q - r o sin fo 
r w = 0 



Vp = (v o sin y 0 + AVj r ) cos fg - (v 0 cos y Q + AV lh )sin fg 
V q = (V 0 sin y Q + AV l)r )sin fg + (v o cos y Q + AV 1>h )cos fg 

V w " ^l,n 

Rjl = cos a) 0 cos D 0 - sin a> 0 sin 0 0 cos i 0 
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APPENDIX A - Continued 


r 12 

= -sin w 0 cos 

- cos u> 0 

sin f2 0 cos i, 

r 13 

= sin J2 0 sin i Q 



r 21 

= cos <x> 0 sin fJ 0 -t 

- sin w 0 cos S7 0 cos iq 

r 22 

= -sin sin £* 0 

+ cos cu 0 

cos S7 0 cos i ( 

r 23 

= -cos 0 0 sin i 0 



r 31 

= sin cu 0 sin i 0 



r 32 

= cos o> 0 sin i Q 



r 33 

= cos i Q 



r x = 

R ll r p + R 12 r q + 

R 13 r w 


r y = 

R 21 r p +R 22 r q + 

R 23 r w 


N 

11 

R 31 r p +R 32 r q + 

R 33 r w 



V x = RllVp + Hi2V q + Ri3V w 

Vy = R21 v p + R 22 v q + R 23 V w 

V z = R 3 lV p + R32 v q + R 33 v w 

T 

r x> r y> r z> V x , Vy, V z , — ► a^, e^, n, 

where the transformation T from Cartesian position and velocity to Keplerian orbital 
elements is by reference 11. 

The geometry of the second maneuver is defined in sketch B. Knowing fj[, 
one can calculate the radius at the maneuver as 

a,(l - e, 2 ) 

If = r 

1 - -1 
1 + cos fj 
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APPENDIX A - Continued 


Z,Z' 


X 


Sketch B 




APPENDIX A - Continued 

eters being considered. If this is not the case, then the assumption is made that the 
sgn (sin = sgn (sin fjj. This assumption is valid except for small values of |sin f°|. 

It is recommended that both signs of sin f® be examined for the smallest AV 2 if this 
is the case. The unit vector ? to the maneuver point can now be computed by 

r x = cos ^o> t + f°) cos - cos ^ sin^cu t + f°)sin 

r y = cos + f°^sin + cos ij. sin^o> t +f°jcos 

r z = sin(cu t + ft) sin i t 

To simplify the further development, rotate to the X r ,Y’,Z f coordinate system (sketch B) 
and compute 


sin 6 - r z 


cos 9 = +(l - sin^ 0 ) 


sin A = 


cos A = 


cos 9 
r x 

cos 9 


1/2 


If cos 0 = 0, there is no solution as will be shown later. Thus, one can write ? as 


r x > = cos 6 


r z < - sin 6 

Now impose the constraint that the latitude of \p equals 5. To this end define a 
unit vector P in the final orbit which is at a true anomaly of yp and is at a latitude 
of 6. Since the true anomaly of ? is fj and the true anomaly of P is \p, the angle 
between ? and P is (f| - Pj- Thus, the minor circle r • P = cos(ff - ipj can inter- 
sect the 5 latitude line at two points (see sketch B) and yield two solutions to the trans- 
fer problem. The transfer with the lowest AV is desired. The vector P can be 
determined from the following system of equations: 
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APPENDIX A - Continued 

r • P = cos (f j - 
P • z = cos (90° - 6) = sin 6 

A A 

p • p = 1 


From the second equation obtain 
P z t = sin 8 

and from the first equation 

P t cos 0 + sin 0 sin 6 = cos 


(fj - *) 


or 


P x' = 


cos^fj - ) - sin 0 sin 8 

cos 8 


If cos 0 = 0, then 0 = ±90° and r is along z. In this case the minor circle defined by 
r • P = cos^f| - t//j is a latitude line and a solution is possible only if cos (f| - ^ = sin 6. 
Since this is very unlikely, this case is not considered and no solution is assumed if 
cos 0 = 0. The third equation yields 



If P^t < 0, then the minor circle r * P = 
line and there is no solution. If 


p y’ > °» 


p y’ = ±U 


2 

F i’) 


cos ^fj - i//J does not intersect the 
then two solutions exist; namely, 


6 latitude 


Both of these solutions must be examined since they both satisfy all the constraints 
placed on the transfer. The solution which has a plane of motion lying closest to the 
transfer orbit costs less and is chosen. To determine this orbit, rotate P in the prime 
system to the original X,Y, Z system; that is 


P x = P x > cos A - Pyt sin A 
P y = P x > sin A + Py» cos A 
P z = sin 6 
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APPENDIX A - Continued 


where Py? > 0 has been arbitrarily chosen. The unit vector Hf normal to the plane 
of the final orbit and in the direction of angular momentum can be found from 

H = Px r , Pxr 


Hf,x" 


p y r z ' P z r y 


l 

P r 

- P r 

sin(fj - i//) 

^z r x 


H f,z 


P r - P r 
Jr x 1 y y x 

— 


The corresponding vector for the transfer orbit is given by 
H t>x = sin sin it 
Ht y = -cos ty- sin ^ 

H t,z = cos i t 

From sketch C it can be seen that the angle 4 which is the angle between the two planes 
of motion at their line of intersection is given by cos 4 = H t • H f , or 

cos 4 = Hf, x + ^y^y + 

At this point the other solution of P is examined, that is Py< < 0, and the value of cos 4 
corresponding to this solution is calculated. The desired P is that one which has the 
largest value of cos 4, or that solution which has a plane of motion lying closest to the 
transfer orbit. Once this value is obtained, the sine of 4 is calculated; that is, 

sin | = x Hf) • r 


or 


sin 


£ = (^yH^z ' n^z r x + ( H t,z H f,x " ^x^z^y + (^x 1 ^ 


,y 


' H t,y H f,x) ] 
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APPENDIX A - Concluded 


The difference in velocities AV 2 between the transfer orbit and the final orbit 
at r remains to be calculated. The geometry of the velocity vectors is shown in 
sketch D, where the following calculations lead to AV 2 : 


Vf = 




sin 


n 


e t sin f° 

(l + 2e t cos f° + e^ ) 

1/2 


1/2 


cos = +(1 - sin 2 y 


0 


V 2 = 




a . n 

r a f ) 




sin y 2 = 


e^ sin f^ 


1/2 


^1 + 2e f cos ff + 

L . 2 X 1 / 2 

cos y 2 = +fl - sin* y 2 j 

AV 2 r = V 2 sin y 2 - Vj_ sin y ^ 

AV 2,h " V 2 cos y 2 cos % - Vj COS 

AV 2,n = V 2 cos y 2 sin * 


AVr 


’ ( AV l 


+ av| h + av|„) 


1/2 
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APPENDIX B 


RELATION BETWEEN PERIOD FACTOR AND VELOCITY FACTOR 


Assume that the elements of two intersecting elliptical orbits are known. The 
impulsive velocity maneuver between them is AV = V2 - V^. Further assume that the 
period factor a required to obtain an intermediate orbit with period Pq, is known; 
that is, 


P a = Pi +a(P 2 - Pj) 


It is desired to obtain a velocity factor a v which when applied to AV yields an inter- 
mediate orbit with period Pq,. The semimajor axis of the at -orbit is 


a a 



The maneuver occurs at the true anomaly f ° and a radius of 


r = 



1 + ej^ cos f° 


Thus, the magnitude of the velocity in the a -orbit is 

r -1I/2 



From the definition of the velocity factor, 

V^ = V^ + « v AV 
Squaring both sides yields 

Va * = (VT + a v AV) • (v^ + a v Av) 


or 


AV 2 a 2 + 2Vj • AVa„ + fvf - V 2 ) = 0 
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APPENDIX B - Concluded 


and the velocity factor is given by the quadratic formula as 

r 2 i 1/2 

■Vj • AV ± (Vj • AV) - AV 2 (vf - V 2 ) 


“v = 


AV* 


where the sign of the radical is so chosen that 


(hi 


+ 1 - a. 


is a 


A similar expression may be obtained for /3 y by replacing the 
eters with similar /3 related parameters. 


minimum. 
a related par am - 


43 



REFERENCES 


1. Gobetz, F. W. ; and Doll, J. R.: A Survey of Impulsive Trajectories. AIAA J., vol. 7, 

no. 5, May 1969, pp. 801-834. 

2. Van Gelder, A.; Beltrami, E.; and Munick, H.: On Minimum-Time Minimum-Fuel 

Rendezvous. J. Soc. Ind. & Appl. Math., vol. 9, no. 3, Sept. 1961, pp. 474-480. 

3. Straly, W. H.: The Phasing Technique in Rendezvous. ARS J., vol. 32, no. 4, Apr. 

1962, pp. 620-626. 

4. Roth H. L.: Use of the Bi-EUiptic Transfer To Accomplish Single Rendezvous. 

Rep. No. TOR-269(4 130-10) -21 (Contract No. AF 04(695)-269), Aerospace Corp., 
Apr. 9, 1964. (Available from DDC as AD 486 040.) 

5. Bender, D. F.: Rendezvous Possibilities With the Impulse of Optimum Two-Impuise 

Transfer. Vol. 16, Pt. One of Advances in the Astronautical Sciences, Norman V. 
Petersen, ed., Western Periodicals Co., Sept. 1963, pp. 271-291. 

6. Eckel, K. G.: Expense Versus Time in Orbital Transfer to Rendezvous. Astrody- 

namics Guidance and Control, P. Contensou, G. N. Duboshin, and W. F. Hilton, eds., 
Gordon & Breach, Inc., 1967, pp. 67-74. 

7. Doll, J. R.: A Method for the Determination of Optimal Multiple -Impulse Rendezvous 

Trajectories. J-970850-7 (Contract NAS 8-24440), United Aircraft Res. Lab., 

Feb. 1970. (Available as NASA CR- 110 19 5.) 

8. McCue, Gary A.: Optimum Two -Impulse Orbital Transfer and Rendezvous Between 

Inclined Elliptical Orbits. AIAA J., vol. 1, no. 8, Aug. 1963, pp. 1865-1872. 

9. Kibler, James F. ; Green, Richard N. ; and Young, George R. : Orbital Trim by Velocity 

Factoring With Applications to the Viking Mission. Proceedings of the National 
Space Meeting of the Institute of Navigation, Inst. Nav., Mar. 1972, pp. 20-27. 

10. Rosenbrock, H. H. : An Automatic Method for Finding the Greatest or Least Value of 

a Function. Comput. J., vol. 3, 1960/1961, pp. 175-184. 

11. McCuskey, S. W. : Introduction to Celestial Mechanics. Addison -Wesley Pub. Co., 

Inc., c. 1963. 


44 



TABTih' 1.- DETERMINATION OF ALINEMENT OB BIT 



















BISECT- FULL TYPE SOLUTION 



+ |l - a v |AV 1 


C = J AP, 

1 to 


C = (I + J)P + KP. + At* + At* w - T 
d o t o t r 

m - 0, 1, 2 , •**, n = 0, 1, 2, *•*, ri^ 


Alinement on 


initial orbit 


- min | |aaC + C^jj 



C = T - At -IP 
3 a ; o o a o 


If e < e. , minimize AV m over n. 


■^(l" C -( ! V g ) C 3- C 4| 


C = J'AF, + At * a - At oa - x + <: 

3 to t o a,t 4^0 


C, = (I + J f )F + At 10 + At oa - T 
4 o o o a.o 


J’ = ci I - I 


Alinement on 


transfer orbit 


/Af + Af - 

/ o o o \ 

l 555® / 


zb -PV £ ) o *' c *ll 


C = J AP 
3 to 


Ci, = (I + J)P + K*P. + At i0 + At* a - c 
4 o t o t a. I 


K' = ci I - I - J 


/Af 10 + Af* a - Af ia \ 
I o t o » 

\ 5So* j 



If e < solution is acceptable. 


C - J AP 
3 to 

C, = (I + J)F + KP . + MP_ + At 10 + At*°+ Ati a - T 
^ o t I o t I a, i 


f io + <° + ifi e - m u 
o t f o 


M = Ci I I - I - J - K 


/ Af 10 + Af *°+ 

I o t 

" \ 3^0 


If e < c^, solution is acceptable. 








TABLE FULL- BISECT TYPE SOLUTION 


Rendezvous 


Alinement on 


initial orbit 


AV_ = AV_ ) + |6 [fiv + |1 - 6„ I flv 


C 1 = L "ft 

C 2 = IP o + (K + L)P t ♦ fit* 0 - T r 

m = 0, 1, 2, 1 ^; n = 0, 1, 2, 


e = min I [mC + C^l 


C_ — T - At -IP 
3 a,o o a o 


If e < £^, minimize AV^ over n. 


Alinement on 


transfer orbit 


e = min ■ | 


C = T - At*° - IF - At* - K»F 
3 a f t o o t t 


/ &f io > 


Af“- Af 
t o 


If £ < E^, minimize AV^ over n. 


Alinement on 


vMi 


CL = At^ a - At? a + L*AP_. - t + T 
3 f t ft & # f a, L 


i<L A 4 - oa 


C, = IP + (K + L' )P. + At + At + At 

4 O t O t t 


Alinement on 


final orbit 


/Af l0 + fif*°* fif 0 *- Af* a \ 

T T IT I ° t t O ] 

a _ 1 ~ K ^ l 360° / 


If £ < solution is acceptable. 


e =; min ( m£ 


C 3 = L iP ft 


-RP) 


c - c, 

3 


Cj = IP + (K + L)F + MP + At*°+ At*°+ At* a - T 
o t f o t fa, 


- \ 5So^ ) 


M = ci I - I - K - L 


If c < £, , solution is acceptable. 

= D 
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TABLE h.- BISECT- BISECT TYPE SOLUTION 




Minimize AV_ over 


A 








TABTJi TRISECT- FULL TYPE SOLUTION 



S 




TABLE 5.’ TRTSKCT- FULL TYPE SOLUTION - Concluded 





TABLE 6.- FULL- TRISECT TYPE SOLUTION 





T AISLE 6.- FULL-TRISECT TYPE SOLUTION - Concluded 
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